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ABSTRACT 

Context. Even the spectroscopic capabilities of today's ground and space-based observatories can not keep up with the enormous flow 
of detections (> 10'' deg"^) unveiled in modern cosmological surveys as: ;) would be required enormous telescope time to perform 
the spectroscopic follow-ups and n) spectra remain unattainable for the fainter detected population. In the past decade, the typical 
accuracy of photometric redshift (photo-z) determination has drastically improved. Nowdays, it has become a perfect complement 
to spectroscopy, closing the gap between photometric surveys and their spectroscopic follow-ups. The photo-z precision for active 
galactic nuclei (AGN) has always lagged behind that for the galaxy population owing to the lack of proper templates and their intrinsic 
variability. 

Aims. Our goal is to characterize the ability of the Advanced Large, Homogeneous Area Medium-Band Redshift Astronomical 
(ALHAMBRA) survey in assigning accurate photo-z's to broad-line AGN (BLAGN) and quasi-stellar objects (QSOs) based on their 
ALHAMBRA very-low-resolution optical-near-infrared spectroscopy. This will serve as a benchmark for any future compilation of 
ALHAMBRA selected QSOs and the basis for the statistical analysis required to derive luminosity functions up to z ~ 5. 
Mettiods. We selected a sample of spectroscopically identified BLAGN and QSOs and used a library of templates (including the SEDs 
of AGN, and both normal and starburst galaxies, as well as stars) to fit the 23 photometric data points provided by ALHAMBRA in 
the optical and near-infrared (20 medium-band optical filters plus the standard JHKs). 

Results. We find that the ALHAMBRA photometry is able to provide an accurate photo-z and spectral classification for ~88% of the 
170 spectroscopically identified BLAGN/QSOs over 2.5 deg^ in different areas of the survey and brighter than m^yg = 23.5 (equivalent 
to rsLOAN ~ 24.0). The derived photo-z accuracy is below 1% and is comparable to the most recent results in other cosmological fields 
that use photometric information over a wider wavelength range. The fraction of outliers (~12%) is mainly caused by the larger 
photometric errors for the faintest sources and the intrinsic variability of the BLAGN/QSO population. A small fraction of outliers 
may have an incorrectly assigned spectroscopic redshift. 

Conclusions. The definition of the ALHAMBRA survey in terms of the number of filters, filter properties, areal coverage, and depth 
is able to provide photometric redshifts for BLAGN/QSOs with a precision similar to any previous survey that makes use of medium- 
band optical photometry. In agreement with previous literature results, our analysis also reveals that, in the < z < 4 redshift interval, 
very accurate photo-z can be obtained without the use of near-infrared (NIR) broadband photometry at the expense of a slight increase 
in the outliers. The importance of NIR data is expected to increase at higher z (z > 4). These results are relevant for the design of 
future optical follow-ups of surveys containing a large fraction of BLAGN, such as many X-ray or radio surveys. 

Key words, cosmology: Observations - galaxies: active - galaxies: distances and redshifts galaxies: evolution - galaxies: high- 
redshift - quasars: general 



* Based on observations collected at the German-Spanish 
Astronomical center, Calar Alto (Almeria, Spain), jointly oper- 
ated by the Max-Planck-Institut fur Astronomie at Heidelberg and the 
Institute de Astrofisica de Andalucia (CSIC). 



1. Introduction 

The role of active galactic nuclei (hereafter AGN) in the for- 
mation of the early structures and their later evolution has been 
reviewed over the past 15 years, becoming a key ingredient of 



1 



I. Matute et al.: Quasi-stellar objects in the ALHAMBRA survey 



galaxy evolution models (e.g. Cattaneo 2002; Menci et al. 2003; 
Croton et al. 2006; Hopkins et al. 2010 and references therein). 
Evidence shows that many, if not all, massive galaxies harbor 
supermassive black holes (SMBHs; e.g. Kormendy & Richstone 
1995). The close interaction between the formation and growth 
of the SMBH and the evolution of its host galaxy were initially 
revealed by: /) the tight correlations between the masses of cen- 
tral SMBHs and the velocity dispersions and luminosities of the 
bulges of many galaxies (Tremaine et al. 2002); //) the remark- 
able similarities in the redshift at which starburst and accre- 
tion activities ocurred and ///) the observation of the so-called 
downsizing effect not only for the galaxy population but also 
for AGN, i.e. the most massive galaxies appear to have assem- 
bled the majority of their stars earlier than lower mass galaxies 
(Cowie et al. 1996; Zheng et al. 2007), while the density of low- 
luminosity AGN peaks at lower z than the more powerful ones 
(e.g. Hasinger et al. 2008 and references therein). Therefore, the 
measure of the space density of AGN with cosmic time not only 
provides information about the relative importance of accretion 
activity to the global energy output in the universe but also places 
important constraints on early structure formation and galaxy 
evolution (e.g. Di Matteo et al. 2005, Hopkins et al. 2010). 

Quasi-stellar objects (QSOs) are the members of the AGN 
family that have particularly high intrinsic luminosities allow- 
ing them to be detected at large distances and to provide unique 
inside into the early history of the AGN-host galaxy interac- 
tion. Moreover, QSOs are potential contributors to the ultraviolet 
(UV) ionizing background (Cowie et al. 2009) and have proba- 
bly played a non-negligible role in the reionization of the uni- 
verse (Fan et al. 2006, Wang et al. 2010). 

The optical selection of QSOs has been performed mainly 
with follow-up spectroscopic observations of color-color se- 
lected candidates (e.g. SDSS, Richards et al. 2002 and 2dF, 
Croom et al. 2004). These observations use slitless or prism 
spectroscopic surveys and poorly efficient flux-limited spectro- 
scopic surveys (e.g. VIMOS-VLT Deep Survey, Gavignaud et 
al. 2006; Bongiorno et al. 2007). A novel technique with re- 
spect to previous selection criteria was introduced by the CADIS 
(Meisenheimeret al. 1998) and COMBO-17 surveys (Wolf et al. 
2003). These photometric surveys used several optical broad- 
and medium-band filters to characterize the nature of the de- 
tected population and derive its photometric redshift (photo-z) 
via the spectral energy distributions (SEDs). The fluxes reached 
by the survey have allowed the study of the high-z QSO popu- 
lation, thus overcoming the problem of QSO incompleteness in 
the redshift interval 2.2 < z < 3.6 (Richards et al. 2002). This 
redshift range is important because it corresponds to the peak 
and the turnover of the observed QSO space density (e.g. Wolf 
et al. 2004). 

Over the past decade, a clearer understanding of the QSO 
evolution has been achieved thanks to a more accurate charac- 
terization of their different SEDs, to a more precise treatment of 
their variability, and to a significant improvement in their photo-z 
determination. This advance has been encouraged by the concep- 
tion of modern cosmological surveys and newly available space- 
based observing facilities (e.g. HST, XMM-Newton, Chandra, 
Spitzer and Herschel among others) that have been able to detect 
a large amount of sources (> 10^). In particular, when a given 
scientific goal does not require detailed knowledge of the spec- 
tral properties of individual objects, properly designed photo- 
metric surveys can provide a highly reliable photo-z and spectral 
classification for each source. These photo-z's are an essential 
complement to the usually small fraction of sources with spec- 
troscopic redshifts in major extragalactic surveys and to more re- 



liably probe the fainter detected population, which is difficult to 
access using current ground-based spectroscopic observatories. 
In addition, photo-z are used to validate uncertain spectroscopic 
redshifts typically obtained for spectra of low signal-to-noise ra- 
tio (S/N) or limited wavelength coverage (e.g. Fernandez-Soto et 
al. 2001). Several computational methods have been developed 
to derive photometric redshifts with increasingly high precision 
(BPZ, HyperZ, LePhare, ZEBRA, AnnZ, EAzY, among others). 
Only recently have photo-z for AGN (Salvato et al. 2009; Luo et 
al. 2010; Cardamone et al. 2010) reached accuracies similar to 
those computed for normal and starburst galaxies (-1-2%; e.g. 
Ilbert et al 2009). A description of the current state of the art 
photo-z computation as well as a detailed performance compari- 
son of various photo-z codes was provided by Hildebrandt et al. 
(2010). 

In this context, we present the analysis of photometric red- 
shift solutions found for a population of spectroscopically iden- 
tified QSOs using the optical and near-infrared (NIR) multi- 
band catalog of the ALHAMBRA survey. The ALHAMBRA 
survey was designed with an optimal filter combination in or- 
der to provide one of the most homogeneous, large, deep, and 
accurate photometric surveys. Given the proposed depth of the 
ALHAMBRA filters of AB~24.5-25, we expect to sample the 
whole QSO LF up to z ~ 4.2 and up to z~6 for sources with 
Mb > -24.2. At the current stage, the survey has mapped 
~2.5deg^ of the sky in seven different regions. The results pre- 
sented here, and the comparison with existing data from other 
cosmological surveys, will prove the capabilities of the survey to 
derive accurate photometric redshifts for the BLAGN/QSO pop- 
ulation. Furthermore, this test will potentially identify redshift 
ranges for which QSO photo-z estimation maybe unreliable or 
QSOs with atypical SEDs that would then be suitable for more 
detailed study. 

This paper is structured as follows. In Sect. 2, we describe the 
current photometric catalog from the ALHAMBRA survey as 
well as the ancillary data available in each of the ALHAMBRA 
fields from other cosmological surveys. This section also intro- 
duces the QSO sample selection. The methodology followed 
during the photo-z determination is discussed in Sect. 3 while 
in Sect. 4 we quantify the precision of our photo-z estimates by 
comparing them with previous results for this type of sources. 
Finally, Sect. 5 discusses the implications of our results and 
planned future analysis. A detailed QSO catalog will be pre- 
sented in a forthcoming paper. Throughout our analysis, we as- 
sume a ACDM cosmology with H,, - 70kms"' Mpc"^, Qa = 
0.7, and Qm = 0.3. Unless otherwise specified, all magnitudes 
are given in the AB system. 

2. Data set 

2.1. Photometric data: The ALHAMBRA survey. 

The ALHAMBRAQ (Advanced, Large, Homogeneous Area, 
Medium-Band Redshift Astronomical) survey provides a 
photometric dataset over 20 contiguous, equal-width, non- 
overlapping, medium-band optical filters (3500 - 9700 A) plus 3 
standard broad-band NIR filters J, H, and Ks over 8 different re- 
gions of the northern hemisphere (Moles et al. 2008). The survey 
aims to understand the evolution of the structures and the differ- 
ent families of extragalactic sources throughout cosmic time by 
sampling a large enough cosmological fraction of the universe. 
This requires precise photometric redshifts for several hundreds 

' http://alhambra.iaa.es:8080 
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Table 1. The ALHAMBRA fields 
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0.015 
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1 


ALH-4 
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+02 12 21 


0.25 


Dec04-May09 


COSMOS 


0.018 


81/77 


3 


ALH-5 


12 35 00.0 


+61 57 00 


0.25 


May05-Jun09 


HDF-N 


0.011 


18/15 


1,4,5 


ALH-6 


14 16 38.0 


+52 25 05 


0.19 


Aug04-Aug09 


EGS-AEGIS 


0.011 


35/33 


1,2,6 


ALH-7 


16 12 10.0 


+54 30 00 


0.50 


Aug04-Jul09 


SWIRE/ELAIS-Nl 


0.007 


11/11 


1,4,7 


ALH-8 


23 45 50.0 


+ 15 34 50 


0.50 


Aug04-Aug09 


SDSS 


0.024 


3/3 


1 


TOTAL QSOs*°> 






2.44 






180/170 







Notes. 

{i)= ALHAMBRA field name; (//,/n')=Central coordinates of tiie field; (/i')= Area covered by each field; (y)= Period between the beginning and the end of the observations in a given field; 
{vO= Name of the cosmological survey for which a paiticular ALHAMBRA field overlaps; (v//)= Mean Galactic reddening along the line of sight derived by Schlegel, Finkbeiner & Davis 
(1998) from the IRAS 100 /jm data; (viii)- Total number of spectroscopically identified QSOs in the field and those within our selection criteria; {tx)= Source of the QSO classification 
and spectroscopic redshift: (1)= Schneider et al. (2010) ; (2)= DEEP webpage, Davis et al. (2003); (3)=Brusa et al. (2010); (4)= Veron-Cetty, M. R & Veron, P (2010); (5)= Barger et al. 
(2008); (6)= C. Willmer prh: comm; (7)= Rowan-Robinson et al. (2008); 

Sum over all the ALHAMBRA fields of the total/selected spectroscopic QSOs. 



of thousands objects and, therefore, a survey with high photo- 
metric accuracy as well as deep and wide spectral coverage over 
a large area. The simulations of Benitez et al. (2009), relating the 
image depth and photo-z accuracy to the number of filters, have 
demonstrated that the filter-set chosen for ALHAMBRA can 
achieve a photo-z precision, for normal and star-forming galax- 
ies, that is three times better than a classical 4-5 optical broad- 
band filter set. The final survey parameters and scientific goals, 
as well as the technical properties of the filter set were described 
by Moles et al. (2008). The survey has collected its data for the 
20+3 optical-NIR filters in the 3.5m telescope, at the Calar Alto 
observator}0, using the wide-field camera LAICA in the optical 
and the OMEGA-2000 camera in the NIR. The full characteriza- 
tion, description, and performance of the ALHAMBRA optical 
photometric system was presented by Aparicio-Villegas et al. 
(2010). The strategy of ALHAMBRA for each run has been to 
observe the fields with the lowest airmasses trying to complete 
the requested integration time for each filter in order to reach 
the planned depth (ABmag~25). In consequence, and depend- 
ing also on the telescope/instrument downtime and weather, the 
time to complete each filter in each field can vary from months 
to several yearfl Therefore, although ALHAMBRA photome- 
try allows us to detect variability, the lack of a common band(s) 
taken during all the observing runs does not allow us to correct 
for its effect. The deep NIR counts in one of the ALHAMBRA 
fields (ALH-8), over ~0.5 deg-, which has a 50% detection effi- 
ciency depth of 7 ~ 22.4, // ~ 21.3, and ~ 20.0 (Vega), have 
been analyzed by Cristobal-Hornillos et al. (2009). Their results 
helped to constrain different type-dependent galaxy evolutionary 
models. 

In this work, we used the seven ALHAMBRA fields for 
which data have been currently observed and reduced. The cen- 
tral coordinates of these 7 fields, the area covered by each, the 
observing epoch, their coincidence with other cosmological sur- 
veys, and their mean Galactic extinction are detailed in Table [T] 
The total number of spectroscopically identified QSOs, as well 
as the fraction selected for our analysis, are given in Col. 7 of 
Table[T] Figure[T]and Table|2]detail the general characteristics of 
the ALHAMBRA 23 filter set. 

The photometric data points used in this work are given by 
the MAG_AUTO measure by SExtractor (Bertin & Arnouts 1996). 



^ http://www.caha.es 

^ Further details of the observations will be provided in a companion 
paper. 




10000 
X (angstroms) 



Fig. 1. Wavelength coverage of the ALHAMBRA filter set for the LAICA camera CCDl 
(thin continuous lines). The position of the most prominent QSO and BLAGN emission 
lines, plus the Lyman limit, are shown evolving with redshift as thick lines. All filters have 
been normalized to unity (see Apaiicio-Villegas et al. 2010 for their true efficiency). Dotted 
lines trace each filter central wavelength (Table(2). 



To avoid the excessive weight of some points in the computa- 
tion of a photo-z, we adopted a minimum photometric error of 
6m - 0.05 (i.e. if a photometric error is smaller than 0.05, it 
is set to 0.05) for the method described in Sect. 3. In agreement 
with the results of other authors (e.g. Bolzonella et al. 2000), 
we found that there is no gain, or even that we obtain poorer 
results for some objects, when we consider 6m < 0.05. The pho- 
tometric data points of each object were corrected for interstellar 
extinction using the values of E(B - V) provided by the maps of 
Schlegel, Finkbeiner & Davis (1998), which are based on IRAS 
lOOjum datfl 

2.2. Sample selection 

Broad-line AGNs (BLAGN hereafter) and QSOs are powerful 
emitters over the entire electromagnetic spectrum. They show 
significant spectral features in the form of intense emission lines 
(with EW ranging from several tens to several thousands of A) 
in the rest-frame UV, optical, and NIR regime. These properties 

http://irsa.ipac.caltech.edu/applications/DUST/ 
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make QSOs easily detectable out to very high redshifts (z ~ 6) 
and perfect candidates to help improve our understanding of the 
accretion mechanisms within SMBHs (Mbh > 10* Mq). They 
also probe the distribution of large-scale structures and the phys- 
ical conditions of the intergalactic medium (IGM). The strong 
features that characterize the QSO optical emission spectrum al- 
low us to test the ALHAMBRA photometry and its ability to 
produce very low resolution spectra. This would provide a cor- 
rect spectral classification and a high-precision redshift estimate 
for the expected population of several thousands of QSOs. 

We selected our initial QSO candidates from a subsample of 
the current ALHAMBRA catalog (v3), which was created using 
the following photometric criteria: 

- A survey quality flag > 0.7. Each ALHAMBRA source was 
flagged with a parameter ("percent-weight" in the catalog) 
that takes into account the total exposure time of a given 
source relative to the maximum for a given field. A low 
value of this flag (< 0.70) indicates that the source is ei- 
ther within a region strongly affected by the dithering pro- 
cess during the observation, contains bad pixels/artifacts, or 
is located near a bright (masked) source. A detailed com- 
parison with the deeper photometry data available for some 
fields from other surveys shows that source detections with 
percent-weight> 0.70 are highly reliable and that a negligi- 
ble fraction of them are spurious. 

- The source must be within the very high confidence mag- 
nitude interval of the survey. The chosen magnitude of ref- 
erence is A678M filter centered on 6789 A and the interval 
is defined by 17.0 < A678M < 23.5. The bright magnitude 
cut ensures that no source saturates any filter, while the faint 
cut avoids sources with photometric errors larger than ^0.2 
magnitudes. 

We decided against the inclusion of a stellarity criteria, as the 
precision of the one derived by the SExtractor (Bertin & Arnouts 
1996) package was valid only for the brighter part of the QSO 
sample (mevs < 22). 

2.3. Spectroscopic data 

To assess the quality and accuracy of the photo-z determi- 
nation for the ALHAMBRA database, we compiled all the 
published or publicly available spectroscopic information for 
BLAGNs/QSOs. The online services and public spectroscopic 
catalogs included the Sloan Digital Sky Survey (SDSSB DR7 
(Schneider et al. 2010), the Deep Extragalactic Evolutionary 
Probe (DEEP/DEEPS Davis et al. 2003), the All-wavelength 
Extended Groth strip International Survey (AEGI^I; Davis et 
al. 2007), the COSMOS0 XMM source catalog (Brusa et al. 
2010), the GOODS-North redshift compilation by Barger et al. 
(2008), the SWIRlfl spectroscopic catalog by Rowan-Robinson 
et al. (2008) and the 13th edition of the Veron-Cetty & Veron 
QSOs catalog (2010; VERONCAlH hereafter). We verified the 
quality of the identifications according to the following crite- 
ria: /) All spectra from the DR7 SDSS catalog were visually 
inspected to determine whether they contained broad-line emis- 
sion; ii) spectra with a high-quality classification flag (flag>3) 

' http://www.sdss.org/ 
* http://deep.ucolick.org 
^ http://aegis.ucolick.org/ 

http://cosmos.astro.caltech.edu/ 
' http://swire.ipac.caltech.edu/swire/swire.html 

http://heasarc.gsfc.nasa.gov/W3Browse/all/veroncat.html 



Table 2. ALHAMBRA filter characteristics 



Name 




FWHM 


AB 


Offset 


{m) 


<o-(m)) 




(jim) 


(jim) 


corr. 




(AB) 




(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


A366M 


0.3661 


0.0279 


0.96 


-0.033 


21.81 


0.08 


A394M 


0.3941 


0.0330 


0.02 


-0.210 


21.60 


0.05 


A425M 


0.4249 


0.0342 


-0.13 


-0.081 


21.65 


0.05 


A457M 


0.4575 


0.0332 


-0.18 


-0.011 


21.62 


0.07 


A491M 


0.4913 


0.0356 


-0.05 


-0.065 


21.54 


0.06 


A522M 


0.5224 


0.0326 


-0.04 


-0.054 


21.48 


0.06 


A551M 


0.5510 


0.0297 


0.01 


0.003 


21.48 


0.07 


A581M 


0.5809 


0.0324 


0.07 


-0.001 


21.39 


0.05 


A613M 


0.6134 


0.0320 


0.13 


0.009 


21.35 


0.06 


A646M 


0.6461 


0.0357 


0.23 


0.006 


21.33 


0.08 


A678M 


0.6781 


0.0314 


0.24 


-0.046 


21.20 


0.06 


A708M 


0.7078 


0.0332 


0.29 


-0.055 


21.14 


0.05 


A739M 


0.7392 


0.0304 


0.34 


0.007 


21.16 


0.06 


A770M 


0.7699 


0.0354 


0.39 


0.000 


21.11 


0.06 


A802M 


0.8020 


0.0312 


0.44 


0.002 


21.06 


0.07 


A829M 


0.8294 


0.0296 


0.48 


0.007 


21.00 


0.08 


A861M 


0.8614 


0.0369 


0.54 


-0.023 


20.91 


0.05 


A892M 


0.8918 


0.0303 


0.50 


0.022 


20.93 


0.08 


A921M 


0.9208 


0.0308 


0.48 


0.028 


20.83 


0.10 


A948M 


0.9482 


0.0319 


0.52 


0.077 


20.65 


0.15 


J 


1.2094 


0.2471 


0.87 


0.104 


20.65 


0.06 


H 


1.6482 


0.2665 


1.38 


0.186 


20.41 


0.08 


Ks 


2.1409 


0.3040 


1.83 


0.155 


20.20 


0.09 



Notes. 

Columns: (1) Filter name; (2) Filter mean wavelength; (3) Filter FWHM; (4) AB-Vega 
magnitude correction: Wab = flivci;a+AB_con"ection; (5) Offsets applied to each filter as 
'^^fmal— '"filicr+offset during the photometric redshift determination (see Sect. 3.3 for details); 
(6) Mean magnitude in each filter band for the spectroscopic sample fSect J2.3l ; (7) Mean 
magnitude errors in each filter band for the spectroscopic sample. 



were selected from the DEEP/DEEP2 and the AEGIS database 
and visually inspected to confirm that they displayed broad-line 
emission; ///) as neither a spectral classification nor a redshift 
quality were given by Barger et al. in the GOODS-N field, valid 
BLAGN/QSO candidates were selected based on their hard X- 
ray luminosity (Lx[2-8 kev]) being brighter that 10"*^ erg s"' ; iv) 
in the COSMOS field, we selected the high-quality public spec- 
tra of BLAGN (flag£3 13, 14, 18, 213, 214, 218; Lilly et al. 
2007), while we considered as bona fide BLAGN/QSO the re- 
maining of XMM sources without public data reported by Brusa 
et al. (2010) as 'bl' based on MMT and IMACS spectroscopy; v) 
all SWIRE sources have high quality spectra (I. Perez-Fournon 
priv. comm.; and vi) all the sources from the VERONCAT were 
considered as bone fide BLAGN/QSOs. No lower redshift or ab- 
solute magnitude cutoff was included in the source selection as 
our goal is to test the efficiency of our method and photome- 
try as good redshift estimators not only for the most powerful 
BLAGNs and QSOs but also for the low redshift Seyfert 1 nu- 
clei, which may provide an important contribution to the total 
light of their host galaxies. In all cases, the match between the 
ALHAMBRA photometry and the spectroscopic catalogs was 
performed using a one arcsec search radius and always identified 
a unique counterpart. The sources of the BLAGN/QSO spectro- 
scopic redshifts for each of the ALHAMBRA fields are: 

- ALH-2: This ~0.5 deg^ field partially overlaps the deep 
strip of SDSS whose DR7 version provides redshifts for 23 
QSOs. The common region of this area with fields of the 

" flags 18 and 218 refer to a spectroscopic redshift computed with a 
single line 
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Fig. 2. Magnitude-red.shift distribution of the selected spectroscopic sample. 

DEEP/DEEP2 survey yields 7 additional sources from their 
data product release 3 (DR;Q) of the DEEP2 spectroscopic 
catalog. Of the 30 spectroscopic redshifts for BLAGN/QSO 
available in the field, 29 (6 from DEEP2 and 23 from SDSS) 
comply with our photometric criteria. 

- ALH-3; This consists of a ~0.25deg^ field area which par- 
tially overlaps with that of the SDSS. The matching between 
ALHAMBRA and SDSS DR7 spectroscopy yields 2 QSO 
redshifts, for which both sources verify our photometric cri- 
teria. 

- ALH-4: This is a -0.25 deg^ field included in the COSMOS 
survey area. The common area contains a total of 81 
BLAGN/QSO redshifts of which 77 (7 from SDSS and 70 
from COSMOS) comply with our photometric criteria. 

- ALH-5: This field covers a ~0.25 deg^ area overlapping that 
of the GOODS-N. There are a total of 18 BLAGN/QSOs 
with spectroscopic redshifts (9 from SDSS, 6 from Barger et 
al. 2008, and 3 from the VERONCAT). Of these, 15 (7 from 
SDSS, 6 from GOODS-N, and 2 from the VERONCAT) 
comply with our photometric criteria. 

- ALH-6: This ~0.19deg2 field is centered on the GROTH 
strip and therefore overlaps with the DEEP2, AEGIS, and 
SDSS surveys. In total, there are 35 spectroscopically iden- 
tified BLAGN/QSO in this field, of which 33 sources (6 
SDSS, 6 DEEP2 and 21 AEGIS) comply with our photo- 
metric criteria. 

- ALH-7: A ~ 0.5 deg ^ field centered on the ELAIS-Nl 
of the SWIRE survey. Sources with spectroscopic redshifts 
and classifications are provided by the catalog of Rowan- 
Robinson et al. (2008), the SDSS spectroscopy and the 
VERONCAT. All 11 BLAGN/QSO found in this field (1 
from SDSS, 8 from SWIRE, and 2 from VERONCAT) com- 
ply with our photometric criteria. 

- ALHAMBRA-8: The 3 QSO specti-oscopic identifications in 
this ~ 0.5 deg^ field are provided exclusively by the QSO 
DR7 catalog of the SDSS. All 3 sources comply with our 
photometric criteria. 

The final spectroscopic catalog of the ALHAMBRA fields 
includes 94% (170/180) of the total numbers of sources spectro- 
scopically identified in the different fields. The ALH4-COSMOS 
field contains ~44% of the sources, followed by the ALH2 and 

http://deep.berkeley.edu/DR3/dr3.primer.htmI 
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Klg. 3- The complete extragalactic (galaxy + AGN) template database used in this work. 
The source and spectral class for each template given in Table [3]and described in Sect. 3.1 . 

ALH6 fields (DEEP2/ AEGIS) which corresponds to -17% and 
19% of the IDs, respectively. Table[T]details the number of iden- 
tified sources in each field, while Figure|2]shows the redshift dis- 
tribution of the selected spectroscopic sample. 

Although the surveys from which the spectroscopic sample 
is extracted encompass a wide range of selection criteria (color 
selection in the SDSS, optical flux-Hmited in zCOSMOS .bright. 
X-ray selected for IMACS and MMT spectroscopy, etc), we 
find their redshift distribution compatible within the errors. 
The mean redshift and Icr dispersion for the sources extracted 
from the different spectroscopic catalogs are: 1.61+0.63 (SDSS), 
1.99+0.68 (zCOSMOSJ^aint), 1.92+0.64 (zCOSMOS_bright), 
1.24+0.45 (MMT), 1.55+0.66 (IMACS), 1.81+0.81 (GOODS- 
N), 1.56+0.88 (DEEP2), 1.56+0.63 (AEGIS), 1.42+0.69 
(SWIRE), and 1.24+0.78 (VeronCat). 

3. QSO photo-z determination 

We used the publicly available template fitting code LePhar^ 
(Arnouts et al. 1999, Ilbert et al. 2006) to estimate redshifts for 
our selected QSOs. The code matches the photometric data of 
each ALHAMBRA QSO source to a library of available tem- 
plates providing the best-fit, spectral classification, and photo- 
metric redshift by means of ax^ minimization process. The min- 

http://www.oamp.fr/people/arnouts/LE_PHARE.html 
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Table 3. Extragalactic template library. 



Index 


SED Name 


Class 


1 


E112 


Elliptical (5 Gyr old) '•"> 


2 


E115 


Elliptical (2 Gyr old) <"> 


3 


Ell 13 


Elliptical (13 Gyr old) 


4 


Arp220 


Starburst 


5 


M82 


" W 


6 


IRAS 20551^250 


"(a) 


7 


IRAS 22491-1808 


"(») 


8 


NGC 6240 


"(») 


9 


SO 


so'"' 


10 


Sa 




11 


Sb 




12 


Sc 


Sc f"' 


13 


Sdm 


Sdm 


14 


Sd 


Sd 


15 


Spi4 


Spiral 


ID 


jcy 1 o 


LjcyierL i .0 


1 T 

I / 


Sey2 


oeyiert 2 ' ' 


18 


IRAS 19254-7245 


Seyfert 2 


19 


QS02 


QSO 2 


20 


hybl.gallO.agn90 


Hybiid 10% SO + 90% QS02 <*' 


21 


hybl_gal20.agn80 




22 


hybl.gal30.agn70 




23 


hybl_gal40_agn60 




24 


hybl_gal50.agn50 




25 


hybl_gal60.agn40 




26 


hybl_gal70_agn30 




27 


hybl_gal80.agn20 




28 


hybl_gal90.agnl0 


Hybrid 90% SO + 10% QS02 


29 


hyb2.gall0.agn90 


Hybryd 10% 122491 + 90% TQSOl <*' 


30 


hyb2_gal20_agn80 




31 


i-.^.i-.o ^riiiA r,r^^nn 

nyDz_gal-5U_agn /U 




JZ 


Vn/KT rrnl/ln irrn/^rt 

iiyDZ_gai'+u_a.gnou 




33 


hvhl anISO aanSO 




34 


hvh? aalfin atin40 




35 


hvh? fial7n apn'^O 




36 


hyb2_gal80_agn20 




37 


hvh? CTnlOn aanin 


Hvhi-vd 90% T27491 + 10% TOSO! 


JO 




OSO Hnw himinncitv^ 
V^l5U l^lUW-lUIIllIUJalLy ^ 






nSO ('hiffh liiminniiitv^ 
l^OW l^Illgll-IUlIUIlUalLy ^ 


40 


TQSOl 




41 






42 


syth cjso — 25 


OSO svnthetic^'^' 


43 


syth_qso-0.50 




44 


syth_qso-0.75 




45 


syth_qso-1.00 




46 


syth_qso-1.25 




47 


QSO.Cristiani 


QSO f"" 


48 


QSO.VVDS 


QSO composite 


49 


QSO.vandenBerk 


QSO composite 


50 


Mrk231 


BALQSO 



Notes. 

Polletta et al. (2007). <*' Salvato et al. (2009). LePhare tem- 
plate database Cristiani & Vio (1990). <''' Gavignaud et al. (2006). 
Vanden Berk et al. (2001). 



3.1. Template selection 

The list of extragalactic templates used in this work are detailed 
in Table[3]and Figure|3] The selection includes SEDs for QSOs, 
Seyferts, starburst, normal galaxies, and stars. To include low lu- 
minosity BLAGN that are partially or completely dominated by 
their host galaxy light, we adopted the hybrid templates (consist- 
ing of a mixture of QSO and host galaxy SEDs) introduced by 
Salvato et al. (2009). The variety of templates is justified by the 
need to test the ability of the ALHAMBRA survey to differen- 
tiate broad-line AGN emission from that of other extragalactic 
sources or stars, and enable us to do a blind search for these 
sources (Matute et al., in preparation). 

In Tab.|3]and Fig.|3] the templates are organized as: 

- Non-active and starburst galaxies: This includes 3 elliptical 
templates of different ages (2, 5, and 13 Gyr; #1-3), the star- 
burst galaxies Arp220, M82, NGC 6240, IRAS 20551, and 
IRAS 224910 (#4-8) and 7 spirals (SO through Sd; #9-15). 
They are all part of the SED library published by Polletta et 
al. (2007). 

- Obscured BLAGN: Includes the Polletta et al. (2007) com- 
posite templates of a Seyfert 1.8 and a Seyfert 2 and the 
Seyfert-2 IRAS 19254 (#16-18). The high-luminosity ob- 
scured sources are represented by a type-2 QSO (QS02) and 
the BALQSO Mrk231 templates (Polletta et al. 2007; #19 & 
#50). We also added the hybrid template library of Salvato et 
al. (2009) defined by 9 different combinations of a SO and a 
QS02 template (#20-28). 

- QSO and hybrid-QSO templates: Here we considered both 
the high and low luminosity SDSS composites (indices #38, 
#39), two templates from Polletta et al. (2007; QSOl and 
TQSOl with indices #40 and #41), the Cristiani & Vio QSO 
SED (#47), the VVDS mean QSO SED (Gavignaud et al. 
2006; #48), and the mean QSO from Vanden Berk et al. 
(2001; #49) also based on SDSS data. Hybrid templates 
include 9 different combinations of the starburst/ULIRG 
IRAS 22491 template and a QSOl template (Salvato et al. 
2009; #29— #37). As the quality and accuracy of the fit im- 
proved in several cases, we completed the list with a set of 
5 synthetic QSO templates (#42^6), covering continuum 
slopes (v"") from a - -0.25 to a = -1.25 (in steps of 0.25) 
below 1 yum and fixed at c = -0.7 above 1 jjm {LePhare tem- 
plate database and references therein). 

Finally, the stellar template database includes 131 spectra 
from the Pickels (1998) stellar library plus 4 spectra of white 
dwarfs from Bohlin, Colina & Finlay (1995) and 19 additional 
templates from the LePhare stellar library. Stellar templates 
were also included because white dwarfs and F/G stars QSOs 
have similar colors to F/G stars in the z=[2-3] redshift interval. 
Thus, our final database contains 204 templates (50 extragalac- 
tic -I- 154 stellar). During the minimization process. Galactic and 
extragalactic templates were used separately. 



imization process accepts the inclusion of user-supplied priors, 
different extinction laws, and the possibility to apply system- 
atic offsets to the different photometric bands in order to achieve 
the best match between the colors of the sample and those pro- 
vided by the template library. The full capabilities and possibili- 
ties of the LePhare minimization code was extensively discussed 
by Ilbert et al. (2006, 2009). Our final selection of templates, 
adopted reddening laws, priors, and systematic offsets are dis- 
cussed in the following sections. 



3.2. Extinction 

The photometry for several of the sources analyzed in this work 
defines a continuum that strongly deviates from single or even 
double power-laws, most probably because of dust obscuration. 
For many of these types of sources, there is an ongoing debate 
about whether (some of) these red QSOs are obscured by ei- 
ther dust or an intrinsically red continuum (Richards et al. 2003; 

May contain an AGN responsible for 20% of the bolometric flux 
(e.g. Veilleux et al. 2009) 
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Fig. 4. Examples of best-fit solutions assuming a SIVIC extinction law for 6 sources covering a wide range of magnitudes (~ 20 < ^(678 < - 23.5) and spectroscopic redshifts {0.7 < 
z < 2.3). Each panel includes the obsei'ved photometry, associated eiTors, and FWHM for each of the 23 ALHAMBRA filter set (blaclt dots, vertical and horizontal error bar respectively). 
Photometric upper limits are indicated by anows. The continuous line shows the best-fit solution, while the open circles give the expected magnitude from the model corrected from 
systematic offsets. Additional info for each source includes: model name, reduced amount of extinction, the normalized probability distribution as a function of z (Pdz), the spectro-2 
(and its source catalog), and the best-fit photo-e solution. The title of each panel is labeled with the source ID in the ALHAMBRA catalog and the measured magnitude in the m^is filter 
(green dot). 



Young, Elvis & Risaliti 2008). Thus, our multiband template 
fit takes into account the possibility of intrinsic dust obscura- 
tion within the source. We adopted the Small Magellanic Cloud 
(SMC) extinction law (Prevot et al. 1994), which has been shown 
to reproduce the observed reddening for mildly obscured QSOs 
at z < 2.2, for which there are no indications of the Galactic fea- 
ture at 2175A (Hopkins et al. 2004; Richards et al. 2003, York et 
al. 2006). Gallerani et al. (2010) appeared to measure some de- 



viation from the SMC extinction law for higher redshift sources 
which is one reason for adopting alternative extinction laws (see 
below). The attenuation due to dust (Ay) is given as a function 
of the color excess E(B-V) as Ay = Rvy.E{B-V). We assumed 

= 3.1 and a color excess in the range [0,1]. 

Furthermore, as our spectroscopic sample includes lower 
redshift Seyfert 1 nuclei, which may have a strong host-galaxy 
contribution, we considered alternative extinction laws as the 
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dust present in different galaxy types follow extinction curves 
that deviate from that of the SMC. These deviations include vari- 
ations of the steepness in the attenuation curve as in the starburst 
extinction law derived by Calzetti et al. (2000), or the presence of 
a broad bump around 2 175 A as found for the Milky Way (MW; 
Seaton et al. 1979; Cardelli et al. 1989) or the Large Magellanic 
Cloud (LMC; Fitzpatrick 1986). Therefore, to reproduce nor- 
mal galaxy and starburst spectra we also considered in the 
minimization solutions based on the LMC, MW, Calzetti's law, 
and Calzetti's law plus the absorption feature around 2175 A. In 
this case, the minimization process takes into account all possi- 
ble SEDs and extinction laws simultaneously, choosing the best 
suited to each source. These additional extinction laws allow us 
to probe their relevance to the accuracy of the results. 

The light attenuation by the inter-galactic medium (IGM) 
was taken into account internally by LePhare following the 
opacity curves, binned into redshift intervals of Az - 0.1, pub- 
lished by Madau(1995). 

3.3. Systematic offsets 

Photometric redshifts depend strongly on the precision of the 
photometry and the capabilities of the template database to re- 
produce the colors of the source population as a function of z. 
If we were to assume that the selected template database is rep- 
resentative of our source population, then for a given filter the 
average deviation between the observed flux and the best-fit pre- 
dicted flux should be zero for normally distributed uncertainties. 
If this is not the case, a zero-point offset must be applied to the 
photometry derived from the template database when there is a 
non-zero average deviation between the observed and predicted 
fluxes. The ALHAMBRA photometric calibration is based on a 
selection of NGSL stars f ollowing the methodology disc ussed 
in Aparicio-Villegas et al. lAparicio-Villegas et al.l (12010 ). We 
note here that we have not modified this criteria. The computed 
offsets instead, "adapt" the templates to provide a better fit the 
observed photometry. Therefore, alternative templates or object 
selection can (and will) lead to different offsets (e.g. Table 1 by 
Ilbert et al. 2006). 

We used the spectroscopic sample described in Sect, 2 and 
the colors in the filters A457M, A646M, and A829M to com- 
pute these systematic offsets using an iterative approach of find 
the best-fit SED for each source, deriving the mean deviations 
for each filter, applying offsets, re-computing best-fit SEDs, etc. 
The iteration process was halted when the variation in be- 
tween iterations drops below 2%. In general, the procedure did 
not require more than 4 iterations to converge. We found that the 
offsets to be applied are small and agree with the typical pho- 
tometric error for each band in the sample. Table |2] reports the 
values of these corrections as well as the typical photometric er- 
ror for each band in columns 5 and 7, respectively. 

3.4. Priors 

The introduction of important a priori information into the red- 
shift probability distribution function (Pdz) based on Bayesian 
probability can in many cases improve the quality of the solu- 
tions by favoring a particular redshift based on known redshifts 
and/or color distributions (e.g. Benftez 2000). Our analysis only 
makes use of a particular luminosity and redshift range prior and 
does not include any redshift distribution or color information of 
known BLAGN/QSO populations. We restricted the permitted 
absolute magnitudes in the A457M (/14575A) filter between -17 
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Fig. 5. Photo-2 efficiency using several extinction laws (the MEL solution). Top: 
Comparison between the best fit photo-z solution and the measured spectro-:: shows a good 
agreement between the both. The continuous line gives the Zphat = Zspec relation while the 
dashed line represent the boundary between good solutions (green dots) and outliers (red, 
indexed dots) and defined as |A2|/(1 +Zspec) > 0. 15. Central: Distribution of A^/(l + e^p^c) as 
a function of the ALHAMBRA magnitude »?678- The mean magnitude eiTor, per magnitude 
interval of lS.m - 1. of the filter A678M is indicated by the continuous lines. The accuracy 
per magnitude interval (Am = 1) is highlighted by a grey shaded area. This accuracy shows 
a small coiTelation with apparent magnitude. Bottom: This panel shows the contribution of 
each magnitude bin (Am = 1) to the outlier population as filled diamonds connected by a 
dashed line. The magnitude intervals considered and the associated eiTors are indicated by 
the large diamonds. The fraction of outliers with respect to the number of sources in the same 
magnitude intervals are given by asterisks connected by a continuous line. The magnitude 
intervals and the associated eiTors are indicated by the large vertical and horizontal lines. 
In both cases, errors are assumed to be Poissonian and were calculated following Gehrels 
(1986). 



and -28. Absolute magnitudes in this filter are consistent with 
the commonly used broad-band standard filter B. This includes 
not only the typical range where most BLAGN/QSOs are found 
(-28 < M457 < -20) but also the range for host-dominated 
BLAGN and normal galaxies (Salvato et al. 2009; Polletta et al. 
2007; Rowan-Robinson et al. 2008). 
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3.5. Photo-z determination summary 

We used the code LePhare to estimate photometric redshifts for 
170 spectroscopically identified BLAGN and QSOs with high 
quality ALHAMBRA photometry. For the minimization pro- 
cess, we considered the following: 

- A database of 204 templates: 154 stellar and 50 extragalactic. 

- Several extinction laws: MW, LMC, SMC, and Calzetti's 
starburst laws with a color-excess range of E{B - 
y)=[0.0,1.0]. 

- Given our template library, we made a correction to the zero- 
point of the filters that show a non-zero average deviation 
between the observed and best-fit predicted magnitudes. 

- A simple luminosity prior of -17 < Ma457m ^ -28. 

- A redshift space interval of < z < 6 binned in redshifts 
intervals of &=0.04. 

FigureH] shows an example of the excellent agreement be- 
tween the data and the fitted template for 6 sources with a 
wide range of magnitudes (20 < m^evsM ^ 23.5), redshifts 
(0.7 <z< 2.3), and intrinsic extinctions (0.0 < E[B -V]< 0.2). 

4. Results and discussion 

4.1. Plioto-z accuracy 

The efficiency of the photo-z determination is quantified by 
comparing the spectroscopic redshifts (hereafter spectro-z) of 
170 sources in our BLAGN/QSO sample. Photometric redshifts 
are generally characterized by both their accuracy and outlier 
fraction. The accuracy is defined as the standard deviation of 
Az/(1 + Zspec), denoted cr^j(i+.^^^^^, where Az = ^spec ^phot 5 
while the outlier fraction (?/) is defined as the fraction of sources 
with catastrophic solutions (i.e. solutions that are inconsistent 
with the measured spectro-z). In our analysis, we assumed that a 
source is an outlier if |Az|/(l +z) > 0.15. This value was selected 
a priori to be compatible with the cutoff of similar studies (e.g. 
Luo et al. 2010, Salvato et al. 2009, Ilbert et al. 2009, Rowan- 
Robinson et al. 2008). An alternative accuracy estimate that has 
been used by several authors (e.g. Ilbert et al. 2006, Brammer et 
al. 2008) is the normalized median absolute deviation (NMAD) 
defined as 

CTNMAD = 1.48 xmedian|^2^f^i2iilM|. 

This parameter can be directly compared to the standard de- 
viation of Az/(1 + Zspec) in the case of normal distributions and 
has the advantage of being less sensitive to outliers. From now 
on, we use o-nmad as our estimate of the photo-z accuracy. 

We now discuss our results based on the number of extinc- 
tion laws considered in the computation, namely either a sin- 
gle (SMC) extinction law (SEL hereafter) or multiple (SMC, 
LMC, Milky-Way, and Calzetti) extinction laws (MEL here- 
after). Table |4] describes the solutions found for the two sets 
of extinction laws considered. In the case of a SMC extinction 
law, we obtained an accuracy of (Tnmad - 0.009 with a fraction 
of outliers of rj ~ 12% (21 out of 170 sources). A comparison 
between the derived photo-z and the spectro-z is shown in the 
top panel of Fig.|5] The narrow scatter present for the good fits 
(green dots) is highlighted by the distribution of the source red- 
shift accuracy that lies in the range |Az|/(l + Zspec) ^ 0.15 (i.e. 
no outlier region) and shown in Fig.|6] This distribution is well- 
represented by a Gaussian with no measurable bias (centered at 
-0.001) and a cr of ~0.006. Identical results were found when 
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Fig. 6. Uncertainty distribution, Az/(l + z), for the 170 BLAGN/QSOs in our sample 
considering the MEL approach. Non-outlier and outliers are represented by filled and open 
histograms, respectively. The continuous line represents the best Gaussian fit to the observed 
distribution of non-outliers. The number of non-outliers (#), the center, and cr of the best-fit 
Gaussian distribution are indicated. 

we considered several extinction laws during the minimization 
process (ctnmad = 0.009, rj ~ 12%) but, as we see in the follow- 
ing paragraphs, the MEL approach is able to more accurately 
reproduce the SED distribution of the BLAGN/QSO population. 

Besides the ability to provide precise photo-z's, our analysis 
allows us to recover the correct SED for most of the sources. 
Figure|2] presents the distribution of the templates for the MEL 
best-fit solutions. We did not find sources with stellar templates 
that have best-fit solutions (i.e. x^^,eiiar < A'^oz-gso) ^"d the major- 
ity of the sources are best-fitted by pure type-1 QSO templates 
or QSO hybrid templates (ULIRG/QSOl template indices #29 
and above; Sect. 3.1). When we did not take into account the 
outlier fraction of the sources, we found that 95.3% of them 
(142/149) have either a QSO or hybrid-QSO best-fit template, 
1.3% (2/149) are fitted by a QS02 or hybrid-QS02 template, 
and 3.4% (5/149) are fitted with a normal or starburst template. 
Five sources have best-fit solution templates compatible with a 
non-active SED. A closer look at the ALHAMBRA photometry, 
the best fit solution, and the observed spectra (when available) 
revealed that: /) One source might be incorrectly classified as 
BLAGN since both the spectra and the ALHAMBRA photom- 
etry point to an early-type galaxy; //) one source belongs to the 
fainter part of the population (mf,78 = 23.33) and the associated 
errors (Am > 0.2) have diluted any possible BLAGN signature in 
the ALHAMBRA photometry; ///) the other 3 sources have op- 
tical spectra compatible with galaxy templates with different de- 
grees of starforming and post-starforming activities, i.e. no signs 
of broad emission and a continuum with a well-defined 4000 A 
break, but with ALHAMBRA photometry showing some degree 
of AGN activity (blue continuum and indications of faint, typ- 
ical BLAGN line-emission namely of Mg ii , C iii] or C iv at 
the spectroscopic redshift of the source). The incorrect solutions 
found for these 3 sources are probably a consequence of the un- 
certainties expected from the method, particularly regarding the 
chosen template database and the absence of any variability cor- 
rection of the photometry. As a test of the degeneracy introduced 
by the chosen template database, we considered an alternative 
database of only QSOs and hybrid QSO/ULIRG templates (see 
indices #29-39 in Table[3]and Fig.|3]l for the 5 sources with a nor- 
mal galaxy or starburst best-fit solution. We found that: /) for 2 
sources we were unable to recover in this case the photo-z, cast- 
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ing some doubts on the spectral classification of the sources or 
an incomplete template database; /;) a correct photo-z was found 
for the other 3 sources where hybrid templates with a weaker 
BLAGN component (10-20%) were selected by the best-fit so- 
lution. This last case highlights the degeneracy introduced by 
the selected template database for certain host and BLAGN lu- 
minosity regimes. For this small fraction of sources (2%; 3/149), 
the incorrect spectral classification, using the method described 
here, will be taken into account in any statistical analysis of the 
BLAGN/QSO population detected in the ALHAMBRA fields to 
be presented in a forthcoming paper (Matute et al., in prepara- 
tion). 

Furthermore, we note that the photo-z determination de- 
scribed here is able to recover the redshift of the sources in the 
interval 2.2 < z < 3.6, which has been traditionally biased 
against the selection of QSOs because of their similar colors 
to F/G stars. Hence, the photometry and the method described 
here could provide an efficient way of both classifing and deriv- 
ing a reliable photometric redshifts for BLAGN/QSO candidates 
pre-selected, for example, by their X-ray flux. The catalog and 
derived luminosity functions for BLAGN/QSO selected purely 
based on the ALHAMBRA photometry will be presented in a 
companion paper (Matute et al. in preparation). 

The dependence of our results on redshift, apparent mag- 
nitude of the source, and the systematic off'sets applied during 
the photo-z computation are described in the following sections. 
These dependences are similar in the two sets of extinction laws 
unless otherwise specified. 



Ellip/Spiral/ Type-2 Hybrid Hybrid Type-1 QSOs 

Starburst QSO SO/QSOa ULIRG/QSOl 



1 11 11 I 11 , 
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1 ' 
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-,, n n ,,,,,, , 




1 






1 1 









10 30 30 40 50 

Template index 



Fig. 7. Spectral energy distribution for the MEL best-fit solutions. Open histogram takes 
into account all sources, while the shaded histograms consider only the sources with good 
photo-; solutions. We find that none of the ALHAMBRA BLAGN/QSOs is well-represented 
by a stellar template. Of all the extragalactic teinplates considered, the majority (95.3%; 
142 out of 149) of the sources with good solutions (i.e. no outliers) have best-fit templates 
compatible with a pure QSO or hybrid QSO/ULIRG template. 

nitude (bottom panel of Fig. |5j, where ~62% of the outliers have 
»J678 > 21.0. Nevertheless, although some outliers might be pro- 
duced by noisier photometry, other factors might also contribute 
to the catastrophic failures (see Sect. 4.3). 



4.1.1. Dependence on redshift 

The accuracy of the photo-z results is rather independent of the 
redshift with the exception of the interval z=[0.9, 1 .4] (light-grey 
square in Fig.|5t for the SMC results). The presence of only one 
prominent line (Mg ii ) in this interval introduces some aliasing 
that depends on the intensity of the line and how well is sam- 
pled by the ALHAMBRA filters. This small degradation of the 
solution occurs when the peak of the Mg ii line falls within two 
ALHAMBRA filters. The distribution of the oudier fraction of 
the sources (Fig.|5] red dots) follows a bimodal behaviour around 
Zspec ~ 1 .4. Below this redshift (Zspec ^1-4), the minimization pro- 
cess tends to overestimate the photo-z solutions, while it under- 
estimates them at higher redshifts (Zspec > 1 .4). This effect can be 
explained by /) the QSO color/redshift degeneracy (i.e. the de- 
gree of similarity between the colors at different redshifts; e.g. 
Richards et al. 2001) and //) a line misidentification (Croom et 
al. 2004). These degeneracies, still present in the ALHAMBRA 
data but to a much lesser extent than for broadband photome- 
try, are highlighted as grey shaded areas and dot-dashed lines 
in Fig.|5]for the color-color and line misidentification degenera- 
cies. Further details of the origin of these degeneracies are given 
in Sect. l4.3l whei'e we explore the nature of the outlier fraction of 
sources. 

4.1.2. Dependency with apparent magnitude 

As highlighted in the central panel of Fig.|5] we do not find any 
dependence of the accuracy on the apparent magnitude of the 
source but only a clear degradation of the solutions is found at 
fainter magnitudes (mevs > 22.0) caused by the slightly noisier 
photometry (Am67g ~ 0.2 at m(,78 = 23 as indicated by the con- 
tinuous line in the central panel of Fig.|5j- On the other hand, the 
outlier fraction shows a moderate correlation with apparent mag- 



4.1.3. Dependence on photometric offsets 

Table |4] reports the improvement in the fit achieved by includ- 
ing the photometric offsets derived from the ALHAMBRA pho- 
tometry and the template dataset. Although the accuracy is only 
marginally higher, we found that the greatest benefit is the re- 
duction in the numbers of outliers (from ~16% to ~12%). The 
same behaviour has been previously observed by other authors 
(e.g. llbert et al. 2009). 

4.2. Extinction distribution 

The distribution of the color excess E(B - V) required by the 
best-fit solutions is shown in Fig. [8] Considering only the SMC 
extinction law (Fig.|8] left), we found that ~ 62% (92 out of 
149) of the non-outlier sources require some dust extinction 
in their best-fit solution. The mean extinction of the sample is 
{E{B - V)} ~ 0.06, with the majority of these sources (~53%) 
having a color excess compatible with the mean Galactic value 
along the line of sight (E(B - V) < 0.05). Only 5 sources re- 
quire extinctions in the range 0.2 < E(B - V) < 0.4. The solu- 
tions obtained using the MEL approach required a slightly larger 
fraction of extincted sources (79%) with a mean color excess of 
{E(B - V)} ~ 0.09. A significant fraction of sources (-41%) 
with small extinctions {E(B - V) < 0.05) was also found in this 
case. We note here that, since several templates in our database 
are already extincted, the derived values of E(B - V) should only 
be considered as lower limits when these template are chosen by 
the best-fit solution. 

The presence of dust-extincted QSOs is no surprise as the 
SDSS survey has established the existence of a non-negligible 
fraction of QSOs with spectral indices (fy oc v") redder than 
= -1 with extinctions as high as E(B -V)~ 0.5 (Gregg et al. 
2002; Richards et al. 2003). The extinction values found in this 
work for BLAGN/QSO are in good agreement with the extinc- 
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Table 4. BLAGN/QSO photometric redshift results 



Medium-bands + JHKs Medium-bands 
SEL<°> MEL*' SEL*"' MEL'*' 





No offsets 


Offsets 


No offsets 


Offsets 


No offsets 


Offsets 


No offsets 


Offsets 




(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


O" NMAD 


0.010 


0.009 


0.011 


0.009 


0.013 


0.013 


0.014 


0.011 


1 


15.9 


12.3 


14.1 


12.3 


16.4 


15.3 


19.4 


15.3 


Extincted fraction 


46.9 


61.7 


61.6 


79.0 


45.8 


55.6 


71.0 


81.3 


{E(B - V)>"" 


0.045 


0.058 


0.085 


0.092 


0.037 


0.049 


0.063 


0.085 


QSOl-Hybrid fraction''"* 


92.3 


92.6 


95.2 


95.3 


92.2 


91.7 


94.9 


95.1 



Notes. 

Results using the SMC extinction law. '''^ Results using several extinction laws {SMC, LMC, MW, and Calzetti). Fraction of sources requiring extinction. '''' Mean color excess applied 
to the extincted sources. Fraction of sources with pure QSOl or hybrid QSOl/IRAS 22491 templates. 



tion interval derived by Richards et al. (2003; EiB - y)=[0.07- 
0.135]) in order to reproduce the spectral indices of red QSOs 
using templates of QSOs with normal colors from the SDSS. A 
detailed analysis of the selection effects introduced by intrinsic 
dust extinction will be carried out for the future catalog release 
of QSOs detected by the ALHAMBRA survey (Matute et al., in 
preparation). 

4.3. Nature of the outliers 

To fully characterize the photo-z estimate for the ALHAMBRA 
QSOs, it is fundamental to understand the reasons for an in- 
correct photo-z determination. In the following, we describe the 
photometry and nature of the ~12% of the sources (21/170) with 
inaccurate photo-z determinations in an attempt to infer the rea- 
sons for the incorrect solutions. The outlier best-fit templates 
(Fig-IZl open fraction of the histogram), extinction values (Fig.|8] 
open fraction of the histogram), and redshift distribution (upper 
panel, Fig.|5]l are similar to those of the remaining population. 
The most probable reasons for the catastrophic failures are enu- 
merated in the following paragraphs. We only discuss the MEL 
fit approach as it offers the best results, but similar conclusions 
can be drawn from the SEL approach. 

1) The less precise photometry of fainter sources. The ap- 
parent dependence of the outlier fraction on magnitude has al- 
ready been discussed in Sect. 4.1.2. This dependence is almost 
expected as the larger photometric errors of these fainter objects 
(continuous line in the central panel of Fig.|5]l broadens the Pdz 
distribution, increasing the number of peaks and directly influ- 
encing the photo-z results. Nevertheless, only 2 outliers (~10% 
of all the outliers; Fig.|5]#'s 4049 & 6024) have large enough 
photometric uncertainties (Am > 0.2) that can explain their in- 
correct photo-z solution. 

2) Emission-line misclassification. Eight outlier solutions 
(38% of all the outliers; Fig.|5]#'s 2022, 4010, 4061, 5015, 5023, 
5016, 6028 and 6035) misclassify either a single or pair of emis- 
sion lines depending on the redshift of the source and on how 
intense and well-sampled the lines are by the ALHAMBRA 
photometry (see Sect. l4.1.TT l. The expected source distribution 
for single line misclassifications (confusion of the C iii] and 
C IV lines with Mg ii and viceversa) is outlined in the top panel 
of Fig.|5]by the dashed lines. The emission-line pairs Mgii - 
HyS and C iii] - Mg ii can be confused in the redshift inter- 
vals z~[0.6-0.9] and z~[1.5-1.9] owing to the limited resolution 
of the ALHAMBRA spectra (i.e. the ability to resolve the line 
given its intensity and position with respect to the ALHAMBRA 
filters). This latter source of confusion is highlighted as the dark 
grey areas in the top panel of Fig|5] For 6 out of the 8 sources, we 



found secondary solutions in agreement with the spectro-z. For 
these sources, additional a priori information would be needed 
in order to favor one solution over another. A simple a posteriori 
condition that weights the solution at each redshift (Pdz), given 
the source magnitude, was derived from known luminosity func- 
tions of optically selected QSOs. We found that this approach 
does not always reduce the fraction of outliers in our sample 
while, if applied to the whole ALHAMBRA survey, may intro- 
duce a bias. 

3) The possibility of an incorrectly assigned spectro-z. The 
ALHAMBRA photometry and the solutions found for 3 sources 
(-14% of all the outhers; Fig.|5]source #'s 4059, 4081 and 7006) 
seem to indicate that an incorrect spectroscopic redshift has been 
assigned. An example of one of these sources is shown in Fig.|9] 
where the ALHAMBRA low-resolution spectra show a well- 
defined power-law continuum with signs of two or more intense 
emission-lines that do not correspond to the assigned spectro-z. 
A low S/N spectra and/or limited wavelength coverage are typi- 
cal causes of an incorrect assignment of spectroscopic redshifts 
(Fernandez-Soto et al. 2001). Two out of three sources are lo- 
cated in the COSMOS field. These sources spectro-z come from 
VIMOS/VLT observations (Brusa et al. 2010) that are not in- 
cluded in the public release of the zCOSMOS database. Some 
degeneracy in the spectro-z is expected given the limited wave- 
length range covered by the VLT/VIMOS observations of A ~ 
[5500, 9500] A. This means that just a single broad emission-line 
will be visible in certain redshift intervals (e.g. confusion may 



#4059 (ALH4-src6) m 



21.92 




10000 
wavelength [A] 

rig. 9. Example of a source with a dubious spectro-z of 1.39. The ALHAMBRA pho- 
tometry and errors, represented by the blact: dots and crosses, provide a photo-z solution 
{continuous line) of 2.53 with peak emission-lines from L^^, C III] , and C IV clearly visi- 
ble at ~4300, ~540() and -6700 A respectively. Lines and symbols as in Fig.[^ The spectra 
of this source is currently unavailable to the public. 
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Fig. 8- Color excess distribution for the sources that required a certain amount of extinction, E(B-V)>0, to be applied to the best-fit template. Left) Color excess distribution using only 
the SMC extinction law. Right) Color excess distribution using various extinction laws (see Sect. 3.2). Histograms as in Fig. IT] 



arise between Mg ii at z ~ [0.8, 1 .5] and Cm at z ~ [2.2, 2.6]), as 
seems the case for these 3 outliers. Although further confirma- 
tion of their true redshifts is required, these sources prove how 
precise photo-z can complement spectro-z making them an in- 
valuable tool in modern observational cosmology. 

4) Intrinsic variability. Seven sources (-33% of all the out- 
Hers; Fig.|5]#'s 4003, 4012, 4039, 4043, 4054, 6033 and 6042) 
have unusual ALHAMBRA SEDs that are incompatible with the 
template database, indicating that variability is probably playing 
an important role. None of the sources display signs of either 
blending or a close companion that might contaminate their pho- 
tometry. Figure[TO]shows one of these sources and illustrates the 
disagreement in the continuum and emission features between 
the ALHAMBRA photometry and the observed spectra. An es- 
timate of the intensity of the variability can be obtained from the 
5 sources within the COSMOS field using the catalog published 
by Salvato et al. (201 1). All 5 sources display flux variability of 
Am > 0.35 and in some cases as large as ISjn ~ 0.80. The impor- 



#6033 (ALH6-srclO) m^^^ = 20.27 



"1 10" 



Modell: qso-0.5_84 
Model2: pl_qSO.dat 
E(B-V) = 0.20 (Calzelti) 
4- 




10000 
wavelength [A] 

Fig. 10. Example of an object with a strange ALHAMBRA SED probably caused by 
intrinsic variability during the time in which photometry was taken. Lines and symbols as 
in Fig.[3 We also indicate with a green line a possible secondary solution with zphot - 0-6. 
As a comparison, we include the observed spectra with a continuous black solid line. Since 
no flux calibration is provided with the observed spectra, we scale the flux accordingly to fit 
the plot scale. The observed spectra has been smoothed using a 20 pixel box. 



tance of the variability correction for the photo-z computation 
of BLAGN and QSOs has been reported in the past by several 
authors (Wolf et al. 2004, Salvato et al. 2009). Unfortunately, 
owing to the observing strategy of ALHAMBRA, it is impossi- 
ble to set the photometry to a common epoch. 

For only one source (#4083) does its incorrect photo-z solu- 
tions seem unclear as it has a fairly bright (mevg = 18.54) power- 
law spectrum with low variability (Am < 0.2 from Salvato et al. 
2011). 

In summary, we found that ~1Q% of the outliers can be ex- 
plained by their less precise photometry, ~38% are caused by 
emission-line misclassification, ~14% shown clear evidence of 
an incorrectly assigned spectro-z, and the remaining ~33% have 
significant intrinsic variability. 

4.4. Performance comparison with other surveys 

Our results are compared to those already published in the lit- 
erature in other cosmological fields. Table|5] summarizes the 
main characteristics of the photo-z accuracy for several relevant 
cosmological surveys, including the wavelength range covered, 
number of photometric bands used, number of sources in the 
spectroscopic sample used to calibrate the photo-z, the depth of 
the photometry, the photometric code used, the number of tem- 
plates, the final precision in terms of the accuracy (cr), and the 
fraction of outliers (//). The table has been divided into one sec- 
tion for the galaxy and one for the AGN results. This highlights 
the accuracy differences between the two populations, as well as 
the significantly larger (by factors of ~ 3-10) spectroscopic 
samples available to the former because AGN represent only 
< 10% of the extragalactic population. We note here that the 
table highlights simply the differences in the wavelength cov- 
erage and filter type used by each survey. The reader must be 
aware of the caveats in this comparison, namely the different 
depths, source populations, methodology, spectroscopic sample 
used, and the computation of the accuracy estimates. 

The benefit of a continuous optical coverage with medium- 
band photometry is clearly illustrated in Tab.|5] where we 
compare broadband surveys such as SWIRE and HDF-N 
with medium-band surveys like MUSYC, COMBO-17, and 
ALHAMBRA itself. The precision of the results are particularly 
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Table 5. Photo-z accuracy comparison in different cosmological fields 



GALAXIES 



Survey 


# Sources'"* 


# Bands"" 


/l-range 


Depth 


Code 


# Temp.'''' 


(T 




Ref.'''' 


COSMOS 


4148 


30(12) 


UV-IRAC 


'Ib ^ 22.5 


LePhare 


31 


0.007 


< 1(0.15) 


(1) 


GOODS-N/HDF-N 


-130 


7(-) 


UBVIJHK 


/814 < 25 


BPZ 


6 


0.06 


-1.0 


(3,4) 


MUSYC (all) 


2551 


32(18) 


Opt-8.0^im 


i? < 26 


EAzY 


semi-analytical 


-0.010 


-5.0(0.10) 


(6) 


SWIRE (all) 


5982 


6-10 (-) 


0.36-4.5 


r < 23.5 


ImpZ-2 


9 


0.035 


-1.0(0.15) 


(7) 


COMBO-17 


813 


17(12) 


0.35-0.93 /jm 


i? < 24 




PEGASE lib. 


-0.07 


-2.0(0.15) 


(8) 


NDWFS-Bootes 


14448 


13(-) 


UV-Opt-NlR-MIR 


i? < 26 


min. 


4LRT^ 


0.041 


5.0 


(10) 


Lockman Hole 


209 


21 (-) 


FUV-Opt-NlR 


Rc < 22.5 


LePhare 


31 


0.034 


10.0(0.15) 


(11) 



AGN 



XMM-COSMOS' 


236 


30(12) 


UV-IRAC 


i'AB =S 22.5 


Le Phare 


30/31 


0.013 


7.2(0.15) 


(2) 


C-COSMOS^ 


236 


30(12) 


UV-IRAC 


'Ib^22.5 


Le Phare 


30 


0.011 


5.1(0.15) 


(2) 


MUSYC (X-ray) 


236 


32(18) 


Opt-8.0^im 


i? < 26 


EAzY 


1« 


0.012 


12.0(0.12) 


(6) 


CDF-S 


446 


35 (18) 


UV-8.0//m 


z<26 


ZEBRA 


265 


0.01-0.06 


1.4-8.1(0.15) 


(5) 


SWIRE (QSO) 


158 


6-10 (-) 


0.36-4.5 fim 


r < 23.5 


ImpZ-2 


9 


0.093 


-33.0 (0.10) 


(7) 


COMBO-17 (QSO) 


52 


17 (12) 


0.35-0.93 /jm 


i? < 24 




1 


0.007 


17.3(0.10) 


(8,9) 


NDWFS-Bootes 


5347 


13(-) 


UV-Opt-NlR-MIR 


i? < 26 


min. 


4LRTS 


0.05-0.18 


5.0 


(10) 


Lockman Hole 


90 


21 (-) 


FUV-Opt-NlR 


< 22.5 


LePhare 


30 


0.069 


18.9 (0.15) 


(11) 


ALHAMBRA 


170 


23 (20) 


Opt-NIR 


»J678 < 23.5* 


LePhare 


50 


0.009 


12.3(0.15) 


(12) 



Notes. 

Number of sources with spectroscopic redshift used for photo-z calibration. 
'''' Maximum number of photometric bands used for photo-z determination. In parenthesis, the number of those filters that have a nanow or medium passband. 
'■^^ Number of extragalactic templates used for photo-z determination. 

idi pgj-cei^t fraction of outliers. Outlier threshold criteria defined as |Az|/(l + Zspec) in parenthesis. 

References: (l)=Ilbert et al. (2009); (2)=Salvato et al. (2011); (3)= Bem'tez et al. (2000); (4) = Coe et al. (2006); (5)=Luo et al. (2010); (6)=Cardamone et al. (2010); (7)=Rowan- 
Robinson et al. (2008); (8)=Wolf et al. (2004); (9)=Wolf et al. (2008); (10)= Assef et al. (2010); (1 1)= Fotopoulou et al. (201 1); (12)= This work. 

We refer here to the sub-sample of point-like/variable sources from the total XMM COSMOS catalog (the QSOV sample). 

Low Resolution Templates. The photometric redshift code allows for an interpolation between all templates (galaxy + AGN). 

This magnitude coiTesponds to roughly the broadband Sloan filter r - 24. 



higher for the BLAGN population, where the accuracy increases 
by a factor of ~10 and the fraction of outliers decreases by a 
factor of 3 or more. 

A more interesting comparison is that between 
ALHAMBRA and the surveys that include a certain number 
of medium optical filters and different wavelength coverages. 
Since the ALHAMBRA survey is a natural successor to the 
COMBO-17 survey (Wolf et al. 2003), it is mandatory to 
compare the performance of its photometric system with 
that of COMBO-17 in the E-CDFS field. We found that, for 
BLAGN/QSO, ALHAMBRA is only marginally worse than the 
latest recaUbration of the COMBO-17 data (Wolf et al. 2008) in 
terms of accuracy (0.009 vs. 0.007). On the other hand, it has 
a smaller fraction of outliers (12.3% vs. 17.3%). The narrower 
medium-band passbands of COMBO-17 (FWHM of -200 A 
vs. ~300A) may be held responsible for the slight accuracy 
advantage, but also introduce a highly redshift-dependent 
selection function owing to the non-continuous optical coverage 
of the filters, as in the case of ALHAMBRA. This continuous, 
non-overlapping coverage of the optical by the ALHAMBRA 
filters, together with the addition of NIR broadband photometry, 
has helped to reduce the number of catastrophic failures. 
Furthermore, our results are obtained without applying the 
variability correction made to COMBO-17 data. 

As part of the Multiwavelength Survey by Yale-Chile 
(MUSYC; Gawiser et al. 2006), photo-z estimates results were 
published by Cardamone et al. (2010) in the E-CDFS using 
the deeper, more extensive, broad-band photometry and also 
including photometry from 18 medium-band filters taken with 
the Subaru telescope. The performance comparison between the 
ALHAMBRA and MUSYC X-ray (AGN-dominated) popula- 
tion yields identical values. This result is unsurprising given that 
the definition of the MUSYC medium-band photometric system 
is very similar to that of ALHAMBRA, i.e. continuous, non- 



overlapping medium-band filters. Nevertheless, we note that, 
given the nature of the E-CDFS observations, the results ob- 
tained by MUSYC extend to much fainter optical magnitudes 
{R ~ 26). Photometric redshifts estimates are also available in 
this field for the sources detected by the Chandra 2 Ms exposure 
in the 436 arcmin^ CDF-S (Luo et al. 2010). Although the high- 
est accuracy and outlier fraction found by the authors is as low 
as 0.01 and 1.4% respectively, they also found that the strong 
dependence of this result on the spectroscopic sample, used to 
correct their SED library, could increase the dispersion and out- 
lier fraction of the result to cti^mad = 0.059 and 77 = 8.1%. A 
larger dispersion would be expected as their analysis also ex- 
tends to a much fainter population (z < 26) as a result of the 
depth of the Chandra observations (~ 2 x 10"'^ ergcm"^ s"' in 
the 0.5-2.0 keV band). The Luo et al. results have been included 
under the AGN section of Table|5j although given the depth of 
the CHANDRA observations, an important fraction of sources is 
expected to be dominated by intrinsically faint AGN, starbursts, 
or even normal galaxies. 

The COSMOS survey (Scoville et al. 2007) provides the 
largest number of BLAGN/QSOs of our ALHAMBRA sam- 
ple thanks to its extensive spectroscopic follow-up (Lilly et al. 
2007; Brusa et al. 2010 and references therein). This spectro- 
scopic follow-up, coupled with the vast multiwavelength infor- 
mation available in the field, has allowed the computation of 
precise photometric redshift for galaxies and AGN (llbert et 
al. 2009; Salvato et al. 2011). Our photo-z determination and 
that of COSMOS share a similar methodology, the use of the 
same photometric code (LePhare) as well as several templates 
in the dataset. A direct comparison is therefore possible between 
their findings and our results for the BLAGN/QSO population 
(the QSOV sample in Salvato et al. 2011), allowing us to fur- 
ther quantify the ALHAMBRA results through: /) the compari- 
son of the photometric system of both surveys, i.e. a COSMOS 
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Fig. 11. Comparison of the photo-z solutions for the 77 common sources in the XMM-COSMOS field. Left) Difl^erences in the photo-2 solutions found by the ALHAMBRA and 
COSMOS photometry as a function of the visual magnitudes in the ALHAMBRA filter A678M. Both surveys agree for ~ 90% of the sources when (zphot-cos^^phot-alh)/1+^spec ^ 0.15. 
where Zphot-cos 'itid Cphot-alh 'ire the photo-2 solutions found in the COSMOS and ALHAMBRA surveys respectively. Right) Comparison of the source-by-source photo-2 accuracy for 
both sui'veys. The threshold limit for outliers. |Az|/(l + Zspec) ^ 0.15, is indicated by vertical and horizontal dashed lines for the COSMOS and ALHAMBRA surveys, respectively. The 
central grey square represent the region where the two surveys agree, while sources in any of the exterior 4 quadrants (A, B, C, and D) have wrong photo-z solution in both COSMOS and 
ALHAMBRA (crosses). Sources within the vertical or horizontal stripes, defined by the dashed lines, and not within the central square have either /) good-COSMOS and poor-ALHAMBRA 
solutions (squares) or //') good-ALHAMBRA and poor-COSMOS solutions (diamonds). 



30 filter-set with a significant number of broad-band filters, in- 
cluding NUV/GALEX and mid-IR/IRAC photometry but also 
12 medium-band filters and 2 narrow -band filters, against the 
23 ALHAMBRA filter-set dominated by optical medium pass- 
bands but with no info in the NUV or mid-IR; //) the impact of 
the variability correction applied to the COSMOS data. 

We found in general good agreement between our photo- 
z and the fraction of BLAGN/QSOs published by Salvato et 
al. (2011) located within the ALHAMBRA field (77 sources in 
~12.5% of the COSMOS area). The results of this comparison 
are presented in Fig.[TT] The left panel shows the overall agree- 
ment between both surveys, with 90% of the sources (69/77) 
lying within (zphot-cos - Zphot-alh)/1 + zspec < 0.15 over the 
entire magnitude range of the sample. However, this comparison 
does not establish whether or when any of the two photomet- 
ric estimates provide an accurate photo-z. This information is 
given in the right panel of Fig.[TTl where we compare the accu- 
racy found by COSMOS and ALHAMBRA source by source. 
The threshold limits for outliers {ctnmad ^ 0.15) are represented 
by vertical and horizontal dashed lines for the COSMOS and 
ALHAMBRA, respectively. The central grey square contains the 
sources for which the two surveys photo-z estimates agree with 
the measured spectro-z (~84%; 65/77). For 8 of the 12 sources 
in disagreement, COSMOS provides a more accurate solution 
{squares in Fig.fTTTi. while 1 source has a more reliable estimate 
in ALHAMBRA {diamond). Of the 8 sources with more accu- 
rate COSMOS solutions, 5 need a strong variability correction 
according to Salvato et al. (2009; see Sect l4.3l l, which may ex- 
plain our wrong estimates, 2 sources have the line misclassifi- 
cation uncertainty described in Sect. l4.3l (#4061 and #4083) and 
another one has large photometric uncertainties in ALHAMBRA 
(#4049). Finally, in both surveys, 3 sources have photo-z esti- 
mates that inconsistent with their spectro-z (quadrants B and C 
in Fig.fTTTi. We note that these photo-z estimates are very similar 
for 3 of the sources (# 4059, 4081 and 4010) and that the red- 
shifts of these sources are based on a single line spectrum. For 
two of them, there seems to be an incorrect measurement of the 
spectro-z (see Sect. 4.3). For the third source, the photo-z solu- 
tion is degenerate, showing two peaks of similar intensity in the 
Pdz distribution, and the method selects the incorrect photo-z. 



To summarize, we observe that the accuracy of the 
ALHAMBRA photometric redshift determination for the 
BLAGN/QSO population is, as expected, much better than any 
broad-band survey and similar to any recent photometric survey 
including medium passband filters. This result is obtained with- 
out the advantage of the larger wavelength coverages of other 
surveys and thanks to the ALHAMBRA filter selection in the 
optical and NIR. The fraction of outliers is comparable to or 
better than other cosmological surveys that include no variabil- 
ity correction. The importance of the variability correction for 
BLAGN/QSOs (Wolf et al. 2004; Salvato et al. 2009) is high- 
lighted by our larger outlier fraction (by a factor of ~2) com- 
pared to a survey such as COSMOS, whose observation strategy 
allowed all its photometry to be set to a common epoch. 

4.5. The impact of the NIR photometry 

The need for NIR photometry in order to increase the accuracy 
of photo-z and reduce the number of catastrophic failures for 
normal and starburst galaxies has been demonstrated by sev- 
eral authors (Connolly et al. 1997; Fernandez-Soto et al. 1999; 
Rowan-Robinson et al. 2003; Ilbert et al. 2006). The subsequent 
improvement (with an increase in accuracy and a reduction in the 
outlier fraction by factors of -2-3; Ilbert et al. 2006) is caused by 
these photometric bands providing a tighter constraint of photo- 
z solutions as they sample the peak emission of the old stellar 
population (~ 1 yum rest-frame) in the z ~[0,1] redshift range. 
At redshifts above z ~ 1.5-2, the peak stellar emission moves 
out the wavelength range covered by the JHKs filters and longer 
wavelength photometry (e.g. Spitzer/IRAC) is required to sam- 
ple this emission. In the case of the AGN population, Salvato 
et al. (2009, 201 1) illustrated that, when only broad-band imag- 
ing is available, sufficiently high photo-z precision and an outlier 
fraction suitable to scientific analysis can only be obtained with 
the addition of NIR photometry. 

We investigated here the importance of NIR photometry to 
the photo-z solutions of BLAGN and QSOs when medium-band 
photometry is available in the optical. To do so, we proceeded 
as described in Sect.[3] without using the J, H, and Ks broad- 
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band filters. The results of the analysis are summarized in Table|4] 
(columns #5 to #8) for the SEL and MEL solutions. We found 
that the impact of the NIR photometry is negligible in terms of 
accuracy (with a marginal decrease of 0.2-0.3%) but does have 
influence on the fraction of outliers, which increases by 3-5% 
depending on the set of extinction laws considered during the 
fit. The larger outlier fraction is dominated by sources with de- 
generate solutions (showing two or more peaks in the Pdz distri- 
bution) that can be constrained by the NIR photometry. On the 
other hand, the distribution of the SEDs and the required extinc- 
tion is also very similar to the solution with and without JHKs 
photometry. This result is important to unveiling and character- 
izing the BLAGN population of future very lai-ge ALHAMBRA- 
like optical photometric surveys such as the JPAI^ survey or to 
planning medium-band optical photometric follow-ups to char- 
acterize the population of the all-sky X-ray survey to be carried 
out by the future eROSITAH mission (Cappellutti et al. 2011), 
as detailed below. 

We conclude that, although the use of NIR photometry 
helps us to alleviate the fraction of outliers, its importance for 
BLAGN/QSO photo-z determination (in the z - [0, 4] redshift 
range) is significantly smaller than for the normal and star- 
burst population when a continuous coverage of the optical with 
medium-band photometry is available. Nevertheless, the impor- 
tance of JHKs photometry obviously increases with redshift and 
in particular at z > 5, where the limited number of optical detec- 
tions due to the continuum depression bluewards of Lya make 
the NIR bands necessary to constrain the continuum slope of 
BLAGN/QSOs and their photo-z solutions. 

4.5.1. A practical case 

In the future, several surveys will provide large samples of 
BLAGN and QSOs. One of these surveys is to be carried out by 
the eROSITA mission. eROSITA will map all the sky more than 
one order of magnitude deeper in X-ray flux {F{2 - lOkeV] ~ 
10"'^ ergcm"^ s"') than the previous ROSAT All Sky Survey. 
The deep part of the survey will cover ~200 deg- down to a flux 
of F[2 - 10 keV] ~ 4 x 10"'"^ erg cm^^ s"' . In total, the eROSITA 
mission is expected to discover more than 10^ new AGN includ- 
ing large fractions of BLAGN and QSOs given the survey flux 
limits. To aid the follow-up of the AGN population uncovered 
by eROSITA and, in order to obtain reasonable photo-z solutions 
(cnmad ~ 0.08 and rj ~ 22%), Salvato et al. (2011) discussed 
the need for at least JHK photometry when only 4-5 broad bands 
are available in the optical, such as those provided by the very 
wide-field surveys LSSlO and Pan-STARSSQ 

Assuming an ALHAMBRA-like survey covering the entire 
sky, we can estimate the number of BLAGN/QSOs that could be 
photometrically identified in eROSITA following the method- 
ology described here. For this purpose, the XMM-COSMOS 
field serves as a benchmark of the accuracy achievable with 
such an approach by cutting the XMM-COSMOS X-ray catalog 
(Brusa et al. 2010) to the fluxes to be reached by the eROSITA 
all-sky and deep surveys as well as the final magnitude cutoff' 
of the ALHAMBRA photometry (r ~ 25). This X-ray selec- 
tion yields a total of 65 (56 with spectro-z) sources within the 
ALHAMBRA-COSMOS region for the deep selection and 24 
(23 with spectro-z) for the shallow all-sky selection. Photo-z so- 
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lutions were computed for these sources using the ALHAMBRA 
photometry and the results are presented in Fig. [12] The photo-z 
accuracies obtained (ctnmad - 0.013, rj = 8.9% for eROSITA- 
deep; ctnmad - 0.012, rj = 4.3% for eROSITA-shallow) are sig- 
nificantly higher than those for broad-band photometry and sim- 
ilar to the results discussed in Sect.|4] The spectroscopic sam- 
ple is dominated by BLAGN/QSOs (82%, 46 out of 56) for 
which the ALHAMBRA photometry is able to recover a cor- 
rect BLAGN/QSO classification (templates from #29 to #49 in 
Fig.© for 78% (36/46) and 90% (17/19) of the deep and shallow 
samples, respectively. The extrapolation of this efficiency from 
the ALH-4/COSMOS area (0.25 deg^) to the eROSITA deep 
(200 deg^) and all-sky extragalactic (20000 deg^) surveys show 
that an ALHAMBRA-like photometric survey can select and 
provide high accuracy photo-z solutions for more than 4x10^ 
and 1x10^ BLAGNs, respectively. 

Unfortunately, the telescope time that would be required 
to carry out a survey such ALHAMBRA, even in the smaller 
deep regions of eROSITA, is too long for any publicly avail- 
able observatory, i.e. the time to execute an ALHAMBRA- 
like survey in 200 deg^ will be ~10 times longer than that of 
ALHAMBRA itself, roughly 40 years. On the other hand, a dedi- 
cated telescope facility optimized for large-area photometric sur- 
veys would do the job in a fraction of the time. This is the case 
of the planned Javalambre PAU (Physics of the Accelerating 
Universe) Astrophysical Survey (JPAS, Benitez et al. 2009). 
This photometric survey, will be carried out at the Javalambre 
observatory with a 2.5m telescope and a 5 deg^ FOV camera, 
and will sample 8000 deg^ of the northern sky over four years, 
starting in late 2013. The JPAS survey will take the photo-z pre- 
cision to the next level with an expected accuracy of o" = 0.003 
thanks to its 54 narrow-band optical filters (FWHM-IOOA). The 
classification and photo-z estimates of BLAGN and QSOs will 
be significantly higher in the common region covered by JPAS 
and eROSITA with respect to ALHAMBRA. This is due to the 
narrower JPAS filter set and its ability to correct for variabil- 
ity. On the basis of the above ALHAMBRA/eROSITA results, 
we expect JPAS to select and provide very accurate photo-z for 
at least 4x10^ BLAGN/QSOs in the JPAS/eROSITA common 
area. 

5. Conclusions and future work 

We have explored the ability of the ALHAMBRA survey pho- 
tometry to assign very accurate photometric redshifts to a pop- 
ulation of BLAGNs and QSOs. We have achieved a precision 
better than 1% using a catalog of 170 spectroscopically iden- 
tified BLAGN/QSO in the ALHAMBRA fields. This preci- 
sion is similar to the previously published photo-z accuracy for 
BLAGN/QSOs in surveys that make use of medium-band opti- 
cal photometry. This has been possible despite the limited wave- 
length coverage of ALHAMBRA relative to other surveys and 
thanks to its photometric system definition (consisting in 20 con- 
tinuous, non-overlapping optical medium-band filters plus the 3 
NIR broadbands J, H, and K,). 

We have used the publicly available code LePhare to de- 
rive our photo-z solutions by means of fitting our photometry 
to a template database of normal and starburst galaxies, type- 
1/2 Seyfert, QSOs, and stars. Our treatment has included a cor- 
rection for systematic offsets between the different photometric 
bands and the template database. Galactic extinction, absorption 
by IGM, and the possibility of intrinsic reddening by adopting 
various extinction laws. In addition to the excellent accuracy of 
out photo-z, our analysis has demonstrated that: 
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Fig. 12. Photometric and spectroscopic redshift comparison for the 
ALHAMBRA/XMM-COSMOS sources with eROSITA deep survey fluxes 
(F[0.5 - 2keV] > 4 X 10"'^ ergcm"- s"'). Symbols color-coded as in Fig|5] Open 
circles indicate the sources with a X-ray flux above F[0.5 - 2 keV] > 4 X 10"'^ ergcm"^ s"' 
(eROSITA all sky). Crosses highlight sources with spectroscopic BLAGN classification, 
while open-squares indicate which sources are classified as such with the ALHAMBRA 
photometry. 



1 . The medium-band filter set used by the ALHAMBRA pho- 
tometry is able to detect the emission and absorption features 
that characterize QSO optical spectra in the redshift interval 
< z < 4. 

2. We have been able to easily differentiate QSO emission from 
stellar emission over the entire QSO redshift range from to 
3. 

3. In the redshift interval < z < 4 and mgvg < 23.5, we have 
characterized the nature of 87.7% of the sources providing 
a correct SED type and robust photo-z. The fraction of out- 
liers (12.3%) shows a clear correlation with magnitude with 
~70% of the outliers being located at mevs > 21.0. 

4. The most probable reasons for our fraction of outliers are: /) 
the faintness or poor photometric quality of some sources, //) 
line misclassification, ///) the intrinsic variability of the AGN 
population (for which the ALHAMBRA photometry applied 
no correction), and iv) the possibility of an incorrect spectro- 
z assignment owing to the limited wavelength coverage of 
the spectra of some objects. 

5. Near-IR photometry is not fundamental to constrain the 
photo-z solutions of BLAGN and QSOs, at least in the red- 
shift interval < z < 4, if the optical regime is cov- 
ered by a continuous medium-band filter set such as that of 
ALHAMBRA. This result is relevant to the design of fu- 
ture optical follow-ups of surveys with a large fraction of 
BLAGN, as in the case of either X-ray or radio surveys. 

We therefore conclude that the analysis performed here val- 
idates the feasibility and accuracy of the photo-z determination 
for a large number of sources with the ALHAMBRA limited (but 
well-defined) set of filters. In the particular case of BLAGN and 
QSOs, our methodology and results suggest the potential abil- 
ity of the ALHAMBRA photometry to detect these sources and 
build a large and unbiased BLAGN/QSO database. The publi- 
cation of such a catalog as well as the precise estimate of its 
inherent incompleteness necessary to derive meaningful statisti- 
cal properties (e.g. number counts and luminosity functions) are 
left to a forthcoming paper 
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